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A  Comparison  of  Data  Reduction  Techniques 
for  Line  Excited  Optoacoustic  Analysis  of  Mixtures 


Scott  B.  Tilden  and  M.  Bonner  Denton 

Department  of  Chemistry 
University  of  Arizona 
Tucson,  Arizona  85721 


ABSTRACT 


The  use  of  discrete  line  laser  excitation  for  the  optoacoustic  analysis 
of  gas  mixtures  holds  promise  as  a  formidable  analytical  tool.  However 
this  technique  presents  challenging  demands  on  the  data  reduction  process  to 
increase  the  reliability  of  the  analytical  results.  The  problems  unique  to 
this  method  which  made  data  reduction  difficult  are  discussed.  Use  of 
simultaneous  equations,  multivariate  least  squares  and  a  new  iterative 
coefficient  weighted  least  squares  approach  are  compared.  Data  indicating 
the  relative  performance  of  these  methods  are  presented. 


Index  headings:  Optoacoustic  Spectroscopic  Mixture  Analysis. 


Introduction 


Optoacoustlc  Spectroscopy  Is  theoretically  a  sensitive  technique  (1) 
capable  of  quantitatively  determining  a  wide  variety  of  gases.  However  there 
has  been  little  work  directed  towards  utilizing  optoacoustlc  spectroscopy  for 
the  analysis  of  gas  mixtures.  This  Is  expecallly  true  In  cases  where  a 
large  amount  of  spectral  overlap  exists  between  many  of  the  components. 

Since  the  ability  to  analyze  gas  mixtures  Is  a  highly  desirable  feature  when 
employing  the  optoacoustlc  technique,  a  comparison  of  various 
data  handling  strategies  is  needed  to  determine  the  limitations  of  these 
strategies  under  actual  analysis  conditions. 


Theory 

The  analysis  of  mixtures  using  any  spectroscopic  technique  becomes 
trivial  when  absorption  regions  of  one  analyte  species  are  not  overlapped 
by  other  analyte  species  In  the  mixture.  The  wide  acceptance  of  the  atomic 
absorption  technique  Is  In  part  due  to  the  fact  that  complicated  mixtures  of 
metals  and  some  anions  can  be  determined  without  any  prior  separation  step. 

The  limited  number  of  situations  where  spectral  overlap  of  atomic  transition  occurs 
can  usually  be  avoided  by  simply  finding  other  atonic  lines  where  spectral  over¬ 
lap  Is  minimized. 

Most  present-day  optoacoustlc  systems  employ  a  low  pressure  IR 
molecular  gas  laser  as  the  light  source.  Due  to  the  nature  of  these  laser 
sources,  only  a  discrete  number  of  laser  lines  may  be  selected.  This 
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can  limit  the  ability  of  the  analyst  to  find  a  laser  line  where  only  one 
component  of  the  gas  mixture  will  absorb,  although  It  Is  possible  to  find 
selected  mixtures  where  spectral  overlap  Is  minimal  (3).  The  spectral 
overlap  problem  Is  further  compounded  In  optoacoustlc  spectroscopy  by  the 
relatively  broad  vibrational-rotational  absorption  bands  when  analyzing  gas 
mixtures  near  atmospheric  pressure. 

The  classical  approach  for  determining  the  component  concentrations  of 
mixtures  where  spectral  overlap  is  significant  is  to  solve  sets  of  simultaneous 
equations: 
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(Eq.  1) 


where  the  general  term  is  the  absorption  coefficient  of  component  i  at 
wavelength  j;  is  the  concentration  of  each  component  i;  and  A^  is  the  measured 
absorption  value  at  each  wavelength  j.  Typically  the  number  of  equations 
equals  the  number  of  unknowns.  Perlmutter,  Shtrikman  and  Slatkine  (4)  used 
the  simultaneous  equations  approach  to  determine  the  presence  of  ethylene 
In  fruit  storage  warehouses.  Other  interfering  gases  Included  C02,  H20 
and  a  wide  range  of  trace  pollutant  gases  which  were  treated  collectively  as 
one  gas  In  the  analysis.  This  study  showed  that  ethylene  could  be  determined 
in  the  presence  of  these  other  gases  under  certain  conditions.  The  authors 
concluded  that  the  accuracy  of  the  analysis  of  ethylene  was  limited  by  the 
accuracy  In  determining  the  spectral  signatures  of  Interfering  gases  (l.e. 
the  absorption  coefficients  at  each  wavelength  used  In  the  analysis).  The 
authors  also  Indicated  that  negative  concentrations  of  C02  were  sometimes 


determined  after  solving  the  four  simultaneous  equations.  A  kinetic  cooling 
argument  was  used  to  explain  this  behavior  (5).  It  can  also  be  shown  that 
negative  results  can  be  explained  by  a  failure  of  the  simultaneous  equations 
approach  to  properly  determine  the  mixture. 

A  possibly  better  mathematical  approach  for  determining  the  constituents 
of  mixtures  Involves  the  use  of  least  squares  multivariate  analysis  (6). 

This  technique  was  first  used  for  mixture  analysis  using  spectroscopic  data 
by  Sternberg,  et  al .  (7)  and  more  recently  by  Warner  et  al .  (8).  Although 
this  technique  is  well  established  it  has  not  been  used  for  optoacoustlc 
mixture  analysis. 

The  least  squares  multivariate  approach  provides  for  the  incorporation 
of  more  equations  than  unknowns  in  the  analysis.  Since  as  much  data  as  is 
available  can  conceivably  be  used,  one  would  intuitively  conclude  that  this 
mathematical  technique  should  be  superior  over  the  simultaneous  equations 
approach.  The  results  of  actual  mixture  analyses  will  confirm  this. 

The  derivation  of  the  multivariate  least  squares  system  has  been  presented 
elsewhere  (7).  The  result  of  the  least  squares  operation  is  the  transformation 
of  a  non-square  set  of  equations  (more  equations  than  unknowns)  into  a  new 
set  of  equations,  where  the  number  of  equations  equals  the  number  of  unknowns. 
The  form  of  this  new  equations  set  is  as  follows: 

C'-C.X  »  C'-A  (Eq.  2) 

where  matrix  C  is  the  original  coefficient  matrix,  C'  Is  the  transpose  of 
C,  X  Is  the  unknown  component  vector  and  A  is  the  measured  or  observed  vector. 
To  determine  X  (the  unknown  concentrations  of  each  component  In  the  mixture) 
Involves  a  matrix  Inversion  step: 

X  ■  [C'-cr’-C'-A 


(Eq.  3) 
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Although  the  matrix  Inversion  step  can  be  tedious,  a  new  matrix,  U,  can  be 
defined: 

U  -  (Eq.  4) 

The  matrix  U  needs  to  be  determined  once  for  each  mixture  system. 

It  will  be  shown  that  the  least  squares  multivariate  analysis  can  give 
more  accurate  results  when  determining  the  component  concentrations  of  a 
mixture.  However  its  most  useful  feature  is  the  concurrent  generation  of 
statistics  which  gives  the  analyst  some  indication  of  the  reliability  of  the 
results.  The  utility  of  the  statistics  generation  step  will  be  apparent 
when  results  of  mixture  analyses  using  the  multivariate  least  squares  approach 
are  presented. 

Additionally  it  is  possible  to  weight  the  least  squares  analysis  to 
achieve  an  even  more  optimum  result.  Weighting  schemes  usually  involve  a 
a  weight  factor  assigned  to  each  observable.  Weighting  factors  can  be 
collected  into  a  weighting  matrix  W.  When  app.ied  to  the  least  squares 
analysis,  these  weighting  factors  will  weight  most  heavily  those  equations  in 
which  the  varlence  of  the  observable  is  a  minimum. 

X  -  [C-W-Cr’-C'-A  (Eq.  5) 

The  optimum  form  of  the  weighting  matrix  W  has  been  studied  (9).  This  study 
showed  that  different  values  of  the  weighting  matrix  W  can  also  heavily  influence 
the  results  of  the  analysis.  It  was  suggested  that  the  optimum  weighting 
matrix  W  should  include  the  reciprocal  variances  of  each  observable. 

Perlumetter  et  al.  (4)  showed  that  the  limiting  factor  in  achieving  a 
reliable  mixture  analysis  with  optoacoustic  spectrscopy  is  the  difficulty  in 
accurately  determining  spectral  signltures  (absorption  coefficients)  of  component 
gases  In  the  mixture.  This  Is  due  primarily  because  direct  overlap  of  a 
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discrete  laser  line  with  the  center  of  the  molecular  absorption  band  is  a 
fortuitous  and  reasonably  rare  event.  More  commonly,  a  laser  line  falls  on  a 
rapidly  decaying  edge  of  the  analyte  absorption  profile.  Very  small  changes 
from  scan  to  scan  in  laser  parameters  can  causes  slight  changes  in  the  lasing 
frequency  and  subsequently  change  the  observed  absorption  coefficient.  Hence 
it  becomes  highly  desirable  to  provide  the  ability  to  weight  the  coefficients 
depending  on  their  statistical  variance  (i.e.  on  how  constant  the  constant  is 
observed  to  be,  experimentally) .  However,  no  strategy  for  directly  weighting 
the  coefficients  used  in  the  least  squares  multivariate  analysis  has  been 
published.  An  iterative  coefficient  weighting  strategy  has  been  devised 
which  accomplishes  this  goal. 

Experimental 


The  iterative  procedure  begins  by  setting  the  weighting  matrix  U  (in 
Eq.  5)  to  ones,  i.e.,  no  weighting  is  initially  used.  The  multivariate  least 
squares  procedure  is  used  to  calculate  an  initial  "guess"  of  the  component  con¬ 
centrations  of  the  mixture.  The  resultant  calculated  component  concentrations 
are  often  in  error.  A  new  weighting  matrix  W  is  calculated  by  determining  a  sum 
of  the  variance  for  each  equation.  The  sum  of  the  variance  is  calculated  by 
summing  the  products  of  each  variance  coefficient  with  Its  calculated  concentration 
value. 

c2jT  “  p2cjlxl  +  °2Cj2X2  +  +  °2CjnXn  ^Eq*  6) 

2 

Where  o  ^  is  the  variance  of  the  coefficient  of  component  1  at  wavelength  j 
2 

and  o  Is  the  sum  of  the  Individual  component  variances.  The  reciprocals  of 

2 

the  calculated  sum  variances  o  jT  are  Introduced  into  matrix  W  and  the  multi¬ 
variate  least  squares  procedure  is  repeated.  A  new  set  of  component  concentrations 
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Is  calculated.  This  Iterative  procedure  Is  repeated  until  the  component  concen- 
tratlons  converge  at  each  successive  Iteration  step. 

It  will  be  shown  that  due  to  the  widely  differing  variances  of  coefficients 
determined  for  many  gases,  a  coefficient  weighting  scheme  Is  appropriate.  The 
utility  of  this  method  will  be  examined  by  comparing  It  against  the  simultaneous 
equations  and  standard  multivariate  least  squares  determinations  of  actual 
gas  mixtures  using  the  optoacoustlc  technique  to  generate  the  mixture  data. 

Tre  ootoacoustic  spectroscopic  system  used  in  this  study  was  comprised  of 
a  C02  lase'-  with  a  cavity  length  of  2.7  meters  and  plasma  tube  lengths  of  2.5 
meters.  T.ne  laser  was  tuned  with  a  75  1/mm  1.5  cm  grating  (Bausch  and  Lomb, 
Rochester,  NY).  A  differential  cell  comprised  of  two  identical  2.5  cm  long 
x  1.25  cr  internal  diameter  cells  connected  by  a  2.5  cm  x  .4  cm  thick  ZnSe  circular 
discs.  Microphones  used  in  this  study  were  Knowles  BW-1789  miniature  microphones 
(0.6  err2  active  area).  A  PAR  125  light  chopper  (PAR-EGG,  Princeton,  NJ)  was 
used  tc  chc?  the  laser  beam  at  a  fixed  frequency  of  260  Hz.  Laser  power  was 
monitored  with  a  KT-2010  pyroelectric  detector  (Laser  Precision  Corp.,  Yorkville, 
NY)  after  passing  the  beam  through  a  diffuser.  The  optoacoustlc  and  pyroelec¬ 
tric  detector  signals  were  monitored  with  PAR  126  and  5101  lock-in  amplifiers 
respectively. 

Gas  mixtures  were  made  by  successive  dilutions  using  gas  syringes  and  gas 
sampling  bulbs.  Pure  gas  samples  of  mulecules  that  are  gases  at  standard 
temperature  and  pressure  were  obtained  by  flushing  gas  bulbs  with  the  pure  gas. 
(Pure  gases  were  obtained  from  Matheson  Gas  Products,  Cucamonga,  CA.)  Gas  samples 
of  molecules  that  are  liquids  at  standard  temperature  and  pressure  were 
obtained  by  allowing  a  small  amount  of  the  liquid  to  equilibrate  In  a  gas-tight 
sampling  bulb.  The  concentrations  of  the  gas  vapor  In  the  bulb  were  calculated 
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by  using  the  Integrated  form  of  the  Clausius-Clapeyron  Equation  and  checking 
with  GC  analysis. 

Optoacoustlc  data  were  taken  by  scanning  the  grating  with  a  stepper  motor 
drive  with  was  under  computer  control.  The  computer  also  supervised  data 
collection  of  the  optoacoustic  and  pyroelectric  laser  power  monitor  signal  from 
the  lock-in  amplifiers.  Normalized  data  (optoacoustlc  signal  divided  by  the 
laser  power  monitor  signal)  at  each  laser  line  was  output  to 
performing  optoacoustlc  scans  of  pure  gas  samples  mixtures.  A  provision 
for  displaying  raw  optoacoustic  scan  data  on  a  Linear  225  strip  chart  recorder 
(Linear  Instruments,  Irvine,  CA)  was  also  provided. 

Results  and  Discussion 


Absorption  coefficients  at  each  laser  line  were  determined  by  scanning  the 
laser  with  a  pure  gas  sample  in  the  differential  optoacoustic  cell.  A  minimum 
of  four  scans  was  used  to  determine  each  coefficient.  Average  coefficients 
at  selected  laser  lines  used  for  a  four-component  mixture  analysis  are 
listed  in  Table  I.  The  standard  deviation  for  each  coefficient  is  also  listed. 
Note  that  the  relative  standard  deviation  for  each  coefficient  of  each  gas  can 
vary  considerably  between  several  different  laser  lines.  Variances  of  the 
relative  standard  deviation  for  each  gas  is  due  to  the  laser  line  overlap 
problem  discussed  previously.  Note  that  generally  the  relative  standard 
deviation  of  coefficients  Is  Independent  of  the  magnitude  of  the  coefficient. 

It  can  be  shown  that  coefficients  with  the  greatest  magnitude  will  have  the 
most  effect  in  determining  the  calculated  unknown  values  when  using  a  mathematical 
technique  based  on  solving  sets  of  linear  equations.  This  supports  the 
contention  that  an  analysis  strategy  based  on  a  coefficient  weighting  scheme 
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might  be  preferred. 

The  first  system  studied  consisted  of  a  furan,  trichloroethylene,  ammonia 
and  freon-12  mixture.  Average  absorption  coefficients  at  the  nine  laser  lines 
used  in  this  study  are  given  In  Table  I.  This  four  component  mixture  was 
treated  initially  using  the  simultaneous  equations  numerical  procedure.  The 
four  optimum  laser  lines  used  In  this  analysis  are  indicated  in  Table  I.  These 
lines  were  chosen  based  on  a  rejection  ratio  calculations  (3).  The  rejection 
ratio  is  defined  as  the  ratio  between  the  absorption  coefficient  of  a  pair  of 
gases  at  a  laser  line.  An  optimum  laser  line  can  be  chosen  for  each  gas  by 
determining  the  line  which  maximizes  the  rejection  ratio  for  this  component 
with  respect  to  other  gas  components  in  the  mixture. 

The  results  using  the  simultaneous  equations  technique  for  a  four  component 
mixture  of  freon-12,  furan,  trichloroethylene  and  ammonia  are  shown  in  Table  II. 
The  calculated  concentration  vatues  for  trichloroethylene,  furan,  and  freon-12 
agree  quite  well  with  the  actual  concentrations.  However  the  calculated 
concentration  of  ammonia  was  negative.  Results  of  the  analysis  of  this  same 
mixture  using  the  multivariate  least  squares  procedure  is  also  shown  in  Table 
II.  This  analysis  procedure  used  data  at  the  nine  lines  listed  in  Table  I. 

Note  that  a  non-negative  concentration  value  for  ammonia  was  determined, 
however  the  standard  deviation  of  the  calculated  concentration  of  ammonia  is 
quite  large  (approximately  100%  RSD) .  Thus  it  could  be  concluded  that  ammonia 
may  be  present  in  this  mixture  but  the  correct  concentration  value  would  be 
in  doubt.  This  illustrates  the  utility  of  the  statistical  generation  aspect 
of  the  least  squares  analysis  procedure. 

The  results  of  the  coefficient  weighted  iterative  least  squares  procedure 
on  the  same  four  component  mixture  are  shown  in  Table  III.  The  results  of  the 
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first  Iteration  are  Identical  with  the  unweighted  least  squares  procedure  since 
the  weighting  matrix  W  Introduced  Into  the  analysis  Initially  contained  ones  In 
all  matrix  elements.  Successive  Iteration  steps  however  tended  to  give  concen- 
tration  values  closer  to  the  actual  concentrations.  By  the  fourth  iteration 
step  all  values  have  converged  and  are  close  to  the  actual  concentrations  of 
the  gas  mixture  in  the  cell.  Ammonia  was  determined  with  only  a  10%  concentra¬ 
tion  error  (30%  RSD)  using  this  analysis  procedure. 

The  results  of  another  four  component  mixture  of  these  same  three  gases 
using  both  the  simultaneous  equations  and  multivariate  analysis  is  shown  in 
Table  IV.  In  this  mixture  the  concentration  of  freon-12  was  increased  by  a 
factor  of  two  and  the  ammonia  concentration  was  decreased  by  a  factor  of  two  from 
the  previous  analysis.  The  simultaneous  equations  and  multivariate  least  squares 
procedure  fail  to  determine  ammonia  as  evidenced  by  the  negative  calculated 
concentration  values. 

This  same  four  compenent  mixture  was  determined  using  the  coefficient 
weighted  least  squares  procedure  as  shown  in  Table  V.  A  positive  value  for 
the  ammonia  concentration  is  evident,  although  the  standard  deviation  Is 
quite  large  (70%  RSD).  The  analyst  could  properly  conclude  that  ammonia  was 
present  in  this  mixture. 

To  further  test  the  effectiveness  of  the  wieghted  least  squares  procedure, 
a  mixture  of  freon-12,  ethylene,  ammonia,  furan  and  methanol  was  added  to  the 
optoacoustlc  cell.  The  analysis  of  this  mixture  could  be  expected  to  be 
difficult  considering  the  magnitude  of  spectral  overlap  between  these  gases,  as 
shown  In  Fig.  1.  The  first  five  component  mixture  tested  consisted  of  90  ppm 
ethylene,  240  ppm  ammonia,  70  ppm  freon-12,  40  ppm  furan  and  30  ppm  methanol. 
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Raw  coefficients  and  standard  deviations  at  each  of  the  20  laser  lines  used  In 
this  analysis  are  shown  in  Table  VI.  Results  using  all  three  techniques  are  shown 
in  Table  VII.  (The  five  laser  lines  used  for  the  simultaneous  equations 
determination  are  indicated  in  Table  VI.)  The  results  of.  the  simultaneous 
equations  analysis  are  poor.  The  concentration  of  freon-12  was  underdetermined 
by  a  factor  of  10  and  the  concentration  of  furan  was  overdetermined  by  100X. 

The  results  of  the  multivariate  least  squares  procedure  are  better,  however 
methanol  was  underdetermined  by  a  factor  of  3  and  freon-12  was  overdetermined 
by  70S. 

Results  using  the  iterative  coefficient  weighted  least  squares  procedure 
(at  the  fourth  iteration  step)  on  the  same  mixture  are  also  shown  in  Table  VII. 
Note  that  all  concentration  values  are  reasonable  a:  the  fourth  iteration  step. 

The  largest  effor  was  in  the  determination  of  the  methanol  and  freon-12  con¬ 
centrations,  which  were  in  error  by  27X  and  40X  respectively. 

A  similar  five-component  mixture  In  which  the  concentrations  of  furan 
and  ammonia  were  decreased  by  a  factor  of  10  and  ethylene  was  Increased  by  a 
factor  of  3  over  the  previous  analysis  was  analyzed.  The  results  using  all 
three  analysis  procedures  are  shown  in  Table  VIII.  The  results  of  the 
simultaneous  equations  procedure  are  very  poor,  for  example  the  concentration 
of  freon-12  was  determined  to  be  negative.  The  results  using  the  multivariate 
least  squares  procedure  are  better.  However  the  concentrations  of  ammonia  was 
overdetermined  by  a  factor  of  5  and  the  concentration  of  methanol  was  underdeter¬ 
mined  by  a  factor  of  2.  The  clear  superiority  of  the  Iterative  coefficient 
weighted  least  squares  approach  is  evident  by  examining  Table  VIII.  All 
calculated  concentration  values  are  reasonably  close  to  the  actual  values  in 


n 


the  cell  by  the  fourth  Iteration. 


Conclusion 


These  investigations  have  demonstrated  the  poor  performance  of  the  simul¬ 
taneous  equations  data  reduction  techniques  as  compared  to  the  multivariate 
least  squares  and  the  iterative  coefficient  weighted  least  squares  approach 
for  determining  gas  mixtures  using  laser  line  excited  optoacoustic  technique. 
The  Iterative  approach  is  observed  to  produce  the  best  results,  particularly 
in  the  most  demanding  situations. 
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TABLE  II.  Comparison  of  Simultaneous  Equations  and 'Multivariate 

Least  Squares  for  the  Analysis  Using  9  Laser  Lines  of  the 
First  Fo'/5-  Component  Mixture. 


Calculated  Concentrations  (ppm) 

Gas  Concentrations 

in  Cell  (ppm) _ Simultaneous  Equations  Multivariate  Least  Squares 


Freon  12  -  270 

394 

329 

±  39 

Trichloro¬ 

ethylene  -  66 

82 

81 

±  6 

Furan  -  475 

403 

301 

±  163 

Ammonia  -  111 

-13 

47 

±  44 

TABLE  III.  Results  at  Each  Step  Using  the  Iterative  Coefficient  Weighted 

ii  ?s*  !Kuares  Analysis  for  the  First  Four  Component  Mixture 
Using  9  Laser  Lines. 


2nd  Iteration 

XI  =  236  ±  50  ppm 
X2  *  71+8  ppm 
X3  =  387  ±  39  ppm 
X4  =  88  +  46  ppm 


4“^'  Iteration 

X’.  =  232  :  46  ppm 
X2  *  72  :  7  ppm 
X3  *  380  :  47  ppm 
X4  =  96  r  32  ppm 


TABLE  IV.  Comparison  of  Simultaneous  Equations  and  Multivariate 
Least  Squares  for  the  Analysis  of  the  Second  Four 
Component  Mixture. 


Calculated  Concentrations  (ppm) 

Gas  Concentrations 

in  Cell  (pom)  Simultaneous  Equations  Multivariate  Least  Squares 


Freon-12 

430 

570 

513 

:  34 

Tri chloroethy- 

lene 

53 

70 

68 

i  5 

Furan 

750 

818 

726 

:  143 

Amr.or.ia 


56 


-64 


-12  s  38 


TABLE  V.  Results  at  Each  Step  Using  the  Iterative  Coefficient  Weighted 
Least  Squares  Analysis  for  the  Second  Four  Component  Mixture. 


Gas  Concentrations  in  Cell:  XT  *  Freon- 12  430  ppm 

X2  *  Trichloroethylene  53  ppm 
X3  =  Furan  750  ppm 
X4  *  Ammonia  56  ppm 


1st  Iteration 

XI  ■  513  ±  34  ppm 
X2  =  68  ±  5  ppm 
X3  *  726  ±  143  ppm 
X4  =  -12  ±  38  ppm 


2nd  Iteration 

XI  =  426  ±  56  ppm 
X2  «  61  ±  5  ppm 

X3  *  787  ±  47  ppm 
X4  =  31  ±  49  ppm 


3rd  Iteration 

XI  =  427  ±  54  ppm 

X2  *  61  ±  5  ppm 

X3  *  784  ±  53  ppm 

X4  »  33  ±  22  ppm 


4th  Iteration 

XI  *  427  +  54  ppm 
X2  *  61  ±  5  ppm 

X3  *  784  ±  53  ppm 
X4  *  33  +  21  ppm 


TABLE  VII.  Comparison  of  Simultaneous  Equations  (SE),  Multivariate  Least 
Squares  (MLS),  and  Iterative  Coefficient  Weighted  Least  Squares 
(ICWLS)a  for  the  Analysis  of  the  First  Five  Component  Mixture 
Using  20  Laser  Lines. 


Gas  Concentrations 


Calculated  Concentrations  (ppm) 


in  Cell  (ppm) 

SE 

MLS 

ICWLS 

Ethylene 

90 

ns 

104  ± 

22 

97  ± 

12 

Ammonia 

240 

229 

228  ± 

84 

235  ± 

44 

Freon-12 

70 

7 

119  ± 

20 

72  ± 

11 

Furan 

40 

34 

30  ± 

20 

34  ± 

8 

Methanol 

30 

34 

11  ± 

10 

29  ± 

8 

aResults  of  the  Iterative  Coefficient  Weighted  Least  Squares  procedure 
taken  at  the  fourth  iteration  step. 


TABLE  VIII.  Comparison  of  Simultaneous  Equations  (SE) ,  .Multi variate 
Least  Squares  (MLS),  and  Iterative  Coefficient  Weighted 
Least  Squares  (ICWLS)3  for  the  Analysis  of  the  Second  Five 
Component  Mixture. 


Calculated  Concentrations  (ppm) 


Gas  Concentrations 
in  Cell  (ppm) 

SE 

MLS 

ICWLS 

Ethylene 

240 

416 

284  ± 

53 

248  ±  14 

Ammonia 

20 

70 

97  ± 

66 

22  ±  6 

Freon-12 

70 

-11 

56  ± 

33 

85  ±  22 

Furan 

4 

28 

4  ± 

3 

3  i  2 

Methanol 

30 

5 

15  ± 

12 

26  ±  5 

aResults  of  the  Iterative  Coefficient  Weighted  Least  Squares  Procedure 
taken  at  the  fourth  iteration  step. 


TABLE  VI.  Raw  Absorption  Coefficients3  and  Standard  Deviations6  at  Each  Laser 

Laser  Line  (100  Transition) 


Gas 

P42C 

P40 

P38 

P36 

P14C 

P10 

R6C 

R8 

R24 

R26 

R28 

PI 

Ethylene 

.06 

.052 

.055 

.06 

.63 

.63 

.06 

.065 

.21 

.17 

.02 

.1 

(.01) 

(.01) 

(.01) 

(.005)  (.1) 

(.1) 

(.002) 

( .002) ( .005) 

(.005) 

(.002)  (.j 

Ammonia 

.66 

.79 

.71 

.67 

.066 

.048 

.93 

.37 

.006 

.006 

.006 

(.04) 

(.064) 

(.021) 

(.02) 

( . 9982 ) ( . 0044 ) ( . 052 ) 

(.12) 

(.002) 

(.002) 

( .002) 

(3 

Freon- 12 

.37 

.37 

.34 

.27 

.016 

.014 

.03 

.015 

.012 

.013 

.014 

(.042) 

(.037) 

(.032) 

(.039) 

( .0049)( .0048) ( .01 ) 

( .005)( .005) 

(.003) 

(.002) 

(.j 

Furan 

.035 

.049 

.056 

.063 

.105 

.105 

.705 

.60 

1.34 

1 .10 

3.76 

i.i 

(.020) 

( .029) 

(.029) 

(.031) 

(.036) 

(.031) 

(.24) 

( .083)( .13) 

(.48) 

(.39) 

<j 

Methanol 

.032 

.027 

.029 

.036 

.050 

.041 

.10 

.19 

.38 

.67 

1.48 

.1 

(.019) 

(.003) 

( . 001 2 ) ( .0066) ( .005) 

(.005) 

( .0025)( . 026 ) ( .022) 

(.030) 

(.07) 

{.i 

d  Raw  Absorption  Coefficients  Measured  with  the  Following  Pure  Gas  Concentrations  In  0pto< 
Freon-12:  700  ppm;  Furan:  400  ppm;  Methanol:  600  ppm. 

b 

Standard  Deviations  Appear  In  Parenthesis  under  each  Coefficient 
c  These  Five  Laser  Lines  were  Used  In  the  Simultaneous  Equations  Analysis  of  the  Five  Com| 


tsa  and  Standard  Deviations6  at  Each  Laser  Line  Used  In  the  Five  Component  Mixture  Study 


,1ne  (100  Transition)  Laser  Line  (020  Transition) 

P10  R6C  R8  R24  R26  R28  P20c  PI  8  PI  6  PI  4  PI  2  P10  R10  R18  R26c 


.63  .06  .065  .21  .17  .02  .04  .06  .02  .01  .04  .005  .005  .005  .005 

(.1)  (.002)  ( . 002 ) ( .005)  (.005)  (.002)  (.002)  ( .003)( .002)( .002)( .002)( .002)( .002)( .002)  (.002) 

.048  .93  .37  .006  .006  .006  .20  .15  .039  .058  .092  .099  .038  .12  .046 

2)(.0044)(.052)  (.12)  (.002)  (.002)  (.002)  (.019)  ( .026)( .017)( .0043)(.021 )( .013)( .004fl ( .022)  (.013) 

1  .014  .03  .015  .012  .013  .014  .02  .017  .014  .017  .018  .026.  .10  .20  .41 

9)( .0048) ( .01 )  ( .005) ( .005)  (.003)  (.002)  (.002)  ( .003)( .002)( .003)( .003)( -004) ( -01S)(  .055)  (.075) 

.105  .705  .60  1.34  1.10  3.76  1.04  1.20  1.17  1.29  2.30  4.54  2.20  1.13  4.40 

)  (.031)  (.24)  ( .003) ( .13)  { .48)  (.39)  (.17)  (.15)  (.14)  (.28)  (.73)  (.93)  (.71)  (.33)  (1./8) 

.041  .10  .19  .38  .67  1.48  .86  .84  .89  .96  1.95  3.57  1.07  .81  .88 

)  (.005)  ( ,0025)( . 026 ) ( .022)  (.030)  (.07)  (.078)  (.11)  (.11)  (.095)(.46)  (.41)  (.18)  (.14)  (.13) 


Ithe  Following  Pure  Gas  Concentrations  In  Optoacoustlc  Cell:  Ethylene:  750  ppm;  Ammonia:  1000  pp:n; 
I;  600  ppm. 


mder  each  Coefficient 

lultaneous  Equations  Analysis  of  the  Five  Component  Mixtures 


